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First principles calculations of the geometric and electronic structures have been 
performed for two coverages (0.25 ML and 0.5 ML) of C on Ni(001) to understand 
the mechanism of the Ni(001) reconstruction induced by carbon adsorption. The 
calculated structural behavior of the system is in a good agreement with experimental 
observations. The calculated path and energetics of the c(2 x 2) - p4g reconstruction 
in Co.5/Ni(001) is provided. A dramatic reduction of the local electronic charge 
on adsorbed carbon is found to occur upon the reconstruction that decreases the 
electron-electron repulsion on C site. This effect together with the formation of 
covalent bonds between C and the second layer Ni atoms, leads to reconstruction of 
Ni(001). 

I. INTRODUCTION 

Nickel is used as a catalyst for many reactions. The presence of carbon in the reaction 
environment may lead to atomic adsorption on the metal surface changing the electronic 
structure and consequently its reactivity 1 . It is thus important to understand the mecha- 
nism by which C adsorbs on Ni surfaces. On the other hand, this system is an excellent 
prototype for addressing fundamental issues such as the effects of the adsorption of light 
elements on the electronic and geometric structures of metal substrate and the character of 
the chemical bonding between the adsorbate and the surface. It is thus not surprising that 
the first observation of C induced reconstruction of Ni(OOl) 2 has attracted much attention 
of experimentalists and theorists who aim to provide increasingly accurate description and 
explanation of the nature of this phenomenoi^ ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 . 

It is well known that when carbon adsorbs on Ni(001) with coverage less than one third 
of a monolayer (ML) it occupies hollow sites and does not change the symmetry of the 
outmost metal layers^. As the coverage exceeds 0.33ML, the surface reconstructs to what has 
become known as "the clock reconstruction"^! Fig. 1 illustrates schematically these two 



2 



geometries, in both of which C atoms occupy the fourfold hollow sites. In the reconstructed 
structure, the topmost Ni atoms are displaced by alternate clock-wise and counter-clock- 
wise rotation around C, forming a geometry of the (2 x 2)p4g symmetry. A large volume 
of additional information has also been obtained from more recent experiments. Scanning 
tunneling microscopy (STM}2. and photoelectron diffraction (PhD) 4 " measurements have been 
performed for non-reconstructed Ni(OOl) with the C coverage less than 0.3ML (C#/Ni(001), 
9 < 0.3ML). Although the quantitative results in the two papers differ, both report a small 
outward radial displacement of the topmost Ni atoms (Nil) around the C adsorbate. On 
the other hand, for 9 > 0.33ML, results of low energy electron diffraction (LEED)^ and 
PhD 4 experiments reveal a large (0.4 - 0.55A) lateral displacement of the Nil atoms and 
penetration of C atoms into the Nil layer that indicate the surface reconstruction. 

The nature of this reconstruction has been a subject of much debate in the past two 
decades. Using a simple force constant model in lattice dynamical calculations. Rahman et 
aljSi have suggested that the reconstruction is driven by phonon softening. Several scenarios 
for the reduction of specific force constants were presented, including those involving coupling 
of C to the Ni atoms in the second layer (N"i2) 2 i?i7i 1 9 which seem plausible since C penetrates 
the surface during the reconstruction. On the other hand, Terborg et al.— argued that the 
reconstruction is caused by C-induced Nil-Nil repulsion rather than by C-Ni2 bonding, a 
conclusion that has been upheld in the recent band structure calculations 14 . Charge transfer 
in the system has also been suggested to be a factor controlling the reconstruction^*- How- 
ever, to our knowledge, the distribution of local charges in the C/Ni(001) system have not 
yet been calculated from first principles. It has also been proposed that surface stress^*^ 
might be responsible for the reconstruction. Ab initio calculations by Alfe et al.— , however 
emphasize chemical bonding effects rather than elastic effects in causing the reconstruc- 
tion. Direct first principles calculations of C induced surface stress on Ni(OOl) 1 ^ also do 
not reveal any specific relationship between surface stress and surface reconstruction. Both 
Alfe et al.— and Hong et al.^ have addressed the issue of geometric and electronic struc- 
tural changes brought about by the adsorption of 0.5 ML of C on Ni(001). Both used the 
plane wave pseudopotential (PWPP) methods with the local density approximation for the 
exchange-correlation potential. The authors of Ref. Q compare the local densities of the 
electronic states (LDOS) for the non-reconstructed and reconstructed structures. However, 

n 

no distinct tractable modification of LDOS is found upon the reconstruction. In Ref. 1131 the 
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valence charge density for the non-reconstructed and reconstructed structures is provided, 
but authors suggest that a more detailed study of the charge density is required. 

The nature of the reconstruction of Ni(OOl) thus remains to be unclear, and its un- 
derstanding requires detailed description of the electronic structure, character of chemical 
bonding and energetics of the system. Since the reconstruction occurs at a certain C cover- 
age 6 ~ 0.33 ML, it is also important to compare and contrast the electronic structures of 
CyNi(OOl) with 6 < 0.33 ML and 6 > 0.33 ML. 

In this paper, we present the results of our calculations of the geometric and electronic 
structure of C/Ni(001) performed for the 0.25 ML and 0.5 ML C coverages within the density 
functional theory with the generalized gradient approximation for the exchange-correlation 
potential. Applying a structural optimization procedure within the PWPP method, we 
reproduce the reconstruction for Co.5/Ni(001) and provide the energetics and path of tran- 
sition from c(2 x 2) to p4g geometry. The valence charge densities and local valence charges 
are calculated for the systems under consideration using more accurate full-potential lin- 
earized augmented plain wave method. Analyzing modification of these characteristics upon 
change in the C coverage and the reconstruction we provide a consistent description of the 
mechanism of the reconstruction. 

II. COMPUTATIONAL DETAILS 

All calculations presented in this paper have been performed within the density functional 
theory with the generalized gradient approximation (GGA) for the exchange-correlation 
functional^- The Ni(001) surface was modelled by a tetragonal supercell consisting of a 7 
layer Ni slab and 11 A of vacuum. For both 0.25 ML and 0.5 ML carbon coverages we have 
used (2x2) two dimensional unit cell which included four Ni atoms per layer. The supercell 
has also contained two C atoms (one atom at each side of the Ni slab) for the 0.25 ML 
coverage and four C atoms (two atoms at each side of the slab) for the 0.5 ML coverage. 
Thus, the supercells have consisted of 30 and 32 atoms for 0.25 ML and 0.5 ML coverages, 
respectively. 

Surface geometry and energetics have initially been calculated using the PWPP method^ 
with ultrasoft pseudopotentials^ used for C and Ni. We set cutoff energies for the plane-wave 
expansion of 300 eV for both clean and C-covered surfaces and used a 5x5x1 Monkhorst- 
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Pack Appoint mesh in the Brillouin zone sampling 18 . The structures were optimized until 
the forces acting on each atom converg ed better than 0.01 eV/A. 

These optimized structures were then used as input for a more detailed analysis of the 
electronic structure, which includes self-consistent calculations of the valence charge density 
and local charges of the systems, using the full-potential linearized augmented plane wave 
(FLAPW) method^ as embodied in the computational code Wien2k22>. Optimized structures 
obtained from PWPP calculations were used as input data for Wien2k which further refined 
the geometries after a few ionic iterations, be consistent, we have used the same supercells 
and k-point mesh in the Brillouin zone for the FLAPW calculations as in the PWPP ones. 
In the FLAPW method, the local charges are calculated through integration over muffin-tin 
(MT) spheres of radius Rmt centered at each atom. To analyze the effect of the adsorbate 
on these specific quantities the set of Rmt should be chosen to be the same for all structures 
under consideration. Ideally, Rmt should be as large as possible without causing the MT 
spheres to overlap. For inner atoms of the Ni slab, including the Ni2 atoms, a choice of 
Rmt = 1.19A is optimal. However, for Nil, Rmt = 1.005A provided more compatibility 
which the shorter C - Ni bond lengths. For C atoms Rmt = 0.773A was needed. To keep 
plane wave cutoff high enough (RK max —7) with the reduced Rmt 's for the surface atoms, 
we used basis sets of about 5500 LAPW 's including 256 local orbitals for the surface with 
0.25 ML coverage and for the reconstructed surface with 0.5 ML coverage and 2800 LAPWs 
including 128 local orbitals for the non-reconstructed surface with 0.5 ML coverage. 

III. RESULTS AND DISCUSSION 
A. Carbon induced structural modification of Ni(OOl) 

1. Effects of 0.25 ML C on Ni(001) geometric structure 

To our knowledge, the case of 0.25 ML C on Ni(001) has not been studied so far using 
first principles calculations based on density functional theory. W3e have carried out such 
calculations of the geometric structure of the system using a structural optimization proce- 
dure that minimizes the forces applied to the atoms. We find our results to be in a good 
agreement with the experimental observations. Namely, the structural optimization does 
not lead to the p4g reconstruction in this system and results in a stable c(2x2) geometry. 
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Even if we initially shift the Nil atoms along the p4g reconstruction path, they revert to 
the c(2x2) geometry during the optimization. We find that the presence of C causes radial 
outward displacement Ar = 0.07A of the Nil atoms from C. This number is in the range of 
the experimental results, which spread from Ar = 0.02 ± 0.03A (for 0.15 ML coverage)^ to 
Ar = 0.15 ± 0.15A(for about 0.2 ML coverage)*. The calculated equilibrium C-Nil separa- 
tion fa = 0.24A is pretty close to d 01 = 0.21 ± 0.07A obtained for C .i 5 /Ni(001) from PhD 
measurements^. 

2. Effects of 0.5 ML C on Ni(001) geometric structure 

We starts our calculation of C .5/Ni(001) with the perfect c(2x2) structure, and calculate 
the total energy (E tot ) of the system as a function of d m . At each step we fix the C position 
and let the rest of atoms relax (move any direction). The result of these calculations is shown 
on the right side of Fig. 2. No reconstruction occurs during structural optimization and a 
E to t minimum is found at doi = 0.53A. Next, we repeat the procedure, but starting with the 
surface Ni atoms (Nil) initially in the p4g reconstructed geometry^. Since the nickel atoms 
are free to relax as the C atoms are lowered, the system finds the stable position of Nil. We 
find that although Nil positions change during optimization depending on doi, the system 
is kept in the p4g structure. 

Since Ni is a ferromagnetic metal, we perform these calculations for both non-magnetic 
and spin-polarized ferromagnetic cases. The values of d m and the lateral Nil displacement 
Axy obtained for the equilibrium (minimum E tot ) structures from these two calculations 
are listed in Tab. |l] together with those measured in experiments. As seen from the table, 
the spin polarization effect on the geometric structure is negligible and both calculated 
structures are in very good agreement with the experimental observations. To understand 
why the magnetic effect on the geometric structure is so small, we calculate local Ni spins in 
the system by integrating spin density over the spheres of 1.2A radius centered at different Ni 
atoms. We find that while the spin on the Ni atoms located at the center of the slab equals 
0.69/jb, it is reduced to 0.39/xb on Ni2 neighboring carbon and diminished to 0.02//^ o n Nil. 
This result indicates that strong covalent C-Ni bonding suppresses spin polarization on the 
Ni atoms surrounding adsorbed C and makes it possible to neglect the magnetic effects on 
the C/Ni(001) surface geometric structure. We thus perform all following calculations for 
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non-magnetic systems. 

The E tot (d i) dependence calculated for non-magnetic p4g structure is plotted on the left 
side of Fig. 2. One can see that the E to t minimum for the p4g structure is about 0.8 eV 
lower than that for the c(2x2) one. At first sight, both dependences shown in Fig. 2 suggest 
that the system has two energy minima: the local one for the c(2x2) structure and the 
global one for the reconstructed p4g structure and there is an energy barrier between them. 
However, our further consideration shows that this is a wrong conclusion. C adsorbed in the 
ideal c(2x2) structure cannot induce forces acting on Ni atoms along any direction, which 
breaks the c(2x2) symmetry. Regular optimization procedures applying the quasi- Newton 
or conjugate gradient methods are thus not able to lead such a system to reconstruction. In 
reality, even if the system appears to be in the c(2x2) structure, atomic vibrations naturally 
break the symmetry that can cause the forces leading to the reconstruction. 

To place these speculations on firmer ground, we perform another calculation starting 
with a system slightly deviated from the (2x2) symmetry. Namely we start with the quasi- 
equilibrium c(2x2) structure obtained above, but with Nil atoms slightly (by O.OlA) shifted 
along the reconstruction path and let all atoms relax. We find that the structural optimiza- 
tion leads Co.5/Ni(001) to the equilibrium p4g structure with parameters (<i i an d Axy) 
practically the same as obtained in the simulation described above. The system thus under- 
goes a c(2x2) - p4g transformation without any energy barrier that is in agreement with 
the conclusion made in Ref~» The path of the Nil and C displacements and the changes 
in E to t obtained during the optimization are shown in Fig. 3. One can see that the Nil 
lateral displacement and C normal shift happen simultaneously, although first Nil move 
faster than C , then, after Axy reaches the value of about 0.4A, it changes slow, while cLqi 
sharply decreases by 0.3A. At the first stage of the reconstruction the total energy of the 
system substantially decreases, while at the second stage E tot (d i) is flat. 

B. Electronic Structure of C/Ni(001) 

Now, when we find that the structural behavior of C/Ni(001) obtained from our calcu- 
lations completely reproduces the experimental observations, we are going to analyze the 
relationship between the geometric and electronic structures of the system to reveal factors 
controlling the clock reconstruction. To this end, using FLAPW method, we calculate the 
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electronic structure for two observed systems: c(2 x 2)C .25/Ni(001) and p4gC .5/Ni(001), 
as well as, for two hypothetical c(2 x 2)C .5/Ni(001) structures: one with quasi-stable geom- 
etry (doi = 0.53A) and another with a reduced C-Nil separation (da = 0.38A). The valence 
charge densities calculated for these structures are plotted along a plane perpendicular to 
the surface in Figs. 4-7. We find strong covalent C-Nil bonds to be formed in all four 
considered structures. These bonds are seen in the figures as "bridges" of electronic density 
between the C and Nil atoms. The c(2 x 2)C . 2 5/Ni(001) and p45fC . 5 /Ni(001) structures 
also have distinct "bridges" between the C and Ni2 atoms, whereas both c(2 x 2)Co.s/Ni(001) 
structures do not. We find thus that the stable structures observed in experiments differ 
from unstable ones by the presence of extra C-Ni2 covalent bonds. The extra covalent bond 
naturally reduces the total energy of the system and increases its stability. 

Now the question emerges why the C-Ni2 bonds are not formed in c(2 x 2)Co.5/Ni(001). 
This bond can occur if the C-Ni2 bond length (lc-Ni.2) is short enough. In Tab. |H]we provide 
values of Ic-Ni2 as well as, the C-Nil bond lengths {Ic-nu) calculated for the four considered 
structures. As seen from the table, in the stable p4gCo.5/Ni(001) and c(2 x 2)Co.25/Ni(001), 
in which distinct C-Ni2 covalent bonds are formed, Ic-Ni2 are much shorter than in the 
unstable structures. However, as we reduce lc-m2 in c(2 x 2)C .5/Ni(001) from 2.33A to 
2. 18 A, the total energy of the system increases (see Fig. 2). To understand the cause of 
such a behavior, we focus on local charges on carbon. The local charges Qmt{C), calculated 
by integrating the valence charge density over the C muffin-tin sphere with Rmt — 0.773 A, 
are listed in Tab. Hill We find that if l C -m2 in c(2 x 2)C . 5 /Ni(001) is reduced by 0.25A, 
Qmt(C) substantially increases. This is caused by reduction of lc-mi that takes place when 
lc~Ni2 decreases: as lc-mi becomes shorter, neighboring Nil atoms compress the electronic 
charge density on carbon. This should enhance the electron-electron repulsion at the C atom 
that may cause the increase in the total energy mentioned above. The reduction of Ic-nu 
happens, because in c(2 x 2)Co.5/Ni(001) there is no room for Nil outward displacement. 
As the clock reconstruction happens, Ic-nu increases and Qmt(C) dramatically decreases. 
It should be noted that Qmt{C) in stable c(2 x 2)Co.2s/Ni(001) is also significantly smaller 
than in the unstable c(2 x 2) structures with the 0.5 ML coverage. However, the C-Nil 
bond length is not the only factor controlling Qmt{C): if the C-Ni2 covalent bond is formed 
some portion of the valence electron density is transferred to the C-Ni2 covalent "bridge". 
This is seen from comparison of Qmt(C) in c(2 x 2)C .25/Ni(001), p4gC .5/Ni(001) and 
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c(2 x 2)C .5/Ni(001) with d m = 0.53A. All three structures have the same Iq-nh = 1.83A 
(see Tab. |TlJ), but different Qmt{C) (see Tab. IIII|) . The C-Ni2 covalent bond is not formed 
in c(2 x 2)C .5/Ni(001), and Qmt(C) in this system is the largest. On the other hand, 
p4gCo.5/Ni(001) has the strongest C-Ni2 covalent bond and the lowest value of Qmt(C). We 
find thus two factors that make the c(2 x 2) - p4g reconstruction favorable in Co.5/Ni(001): 
formation of C-Ni2 covalent bond that reduces the energy of the system and reduction of 
Qmt{C) that prevents strong electron-electron repulsion on carbon. 

Since interaction of dipoles formed in the surface upon adsorption could also affect surface 
geometry^, we calculate the work functions $ of the systems which reflect change in the 
dipole moments. We find, however, that $ correlates rather with the C-Nil separation than 
with in-plane displacements of the Nil atoms. 

In summary, we have performed first principles calculations of the geometric and elec- 
tronic structures of C/Ni(001) with the C coverages of 0.25 and 0.5 ML. The structural 
behavior of the system obtained from the calculations is in a good agreement with experi- 
mental observations. No energy barrier is found on the reconstruction path. We find that 
formation of the C-Ni2 covalent bonds is crucial for stability of C/Ni(001). These bonds 
are present in the stable c(2 x 2)Co.2s/Ni(001) and p4gCo.5/Ni(001) structures and absent in 
unstable c(2 x 2)C .5/Ni(001). Apart from formation C-Ni2 covalent bond, the c(2 x 2) - p4g 
reconstruction in Co.s/Ni(001) is accompanied by a dramatic reduction of the local electronic 
charge on adsorbed carbon that is supposed to decrease the electron-electron repulsion on 
C site. 
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IV. FIGURE CAPTIONS 



Fig. 1 (Color online). Illustration of the c(2 x 2)C .25/Ni(001) (upper panel) and 
p4gC .5/Ni(001) (lower panel) geometric structures. 
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TABLE I: Comparison of structural parameters of Co.5-p4g/Ni(001) calculated in this work with 
those obtained from experiments-^. NM and FM denote non-magnetic and ferromagnetic calcula- 
tions, respectively 

doi(A) Axy(A) 
NM 0.13 0.51 

FM 0.14 0.51 

Exper. 0.1-0.12 0.41 - 0.55 



TABLE II: C-Ni bond lengths calculated for the considered C/Ni(001) structures. 



Ic-mi (A) h-Ni2 (A) 

pig-C . 5 1.83 1.99 

c(2 x 2)-C .25 1-83 2.04 

c(2 x 2)-C .5,d i = 0.38 1.79 2.18 

c(2 x 2)-C .5,d i = 0.53 1.83 2.33 



Fig. 2. Dependence of the total energy versus the C-Nil separation calculated for 
Co.5/Ni(001) in c(2 x 2) and p4g geometries. Zero energy corresponds to C desorption. 

Fig. 3. Reconstruction path and energetics calculated for Co.5/Ni(001). 

Fig. 4 (Color online). Valence charge density calculated for c(2 x 2)C .25/Ni(001) and 
plotted along the plane perpendicular to the surface. 



TABLE III: C local charges and work functions calculated for the considered C/Ni(001) structures. 



Qmt(C) $(eV) 

c(2 x 2-C . 5 ,doi = 0.53A 4.58 5.67 

c(2 x 2)-Co.B,doi = 0.38A 4.67 5.55 

p4g-C .5,doi =0.13A 4.30 5.39 

c(2 x 2)-C .25 4.41 5.35 
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Fig. 5 (Color online). The same as in Fig. 4, but for p4#C . 5 /Ni(001). 

Fig. 6 (Color online). The same as in Fig. 4, but for c(2 x 2)C . 5 /Ni(001) with 
doi = 0.53A. 

Fig. 7 (Color online). The same as in Fig. 4, but for c(2 x 2)C . 5 /Ni(001) with d 0l = 
0.38A. 
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